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Abstract. We have measured the proper motion of the Local Group galaxy IC 10 with the Very Long Baseline Array by 
measuring the position of an ILO maser relative to two background quasars over 4.3 years. After correction for the rotation of 
the Milky Way and IC 10, we obtain a motion of -39+9 pas yr~' toward the East and 31+8 pas yr~' toward the North. This 
corresponds to a total space velocity of 215 ± 42 km s -1 relative to the Milky Way for an assumed distance of 660+66 kpc. 
Assuming that IC 10 and M33, for which also a proper motion measurement exists, are bound to M31, we calculate a lower 
limit for the mass of M31 of 7.5 xl0 u M o . 
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1. Introduction 

1.1. Proper Motions in the Local Group 

Proper motion measurements of Local Group galaxies 
■ are important for our understanding of the dynamics 
and evolution of the Local Group. Presently, measure- 
ments of extragalactic proper motions by optical tele- 
scopes are limited to the most nearby companions of the 



Milky Way, i.e. the LMC 


[Tones, Klemola, & Lin 1994| 


IKallivayalil et al. 2006| 


Pedreros, Costa, & Mendez 2006 1, 


the SMC (|Kallivaya 


il, van der Marel, & Alcock 2006 1, 



the Sculptor dwarf spheroidal galaxy (dShp) 



dSchweitzer et al. 19951 IPiatek et al. 2006>, the Canis 
Major dwarf galaxy ( Dine scu et al. 2005b> . the Ursa 
Minor dSph JPiatek et al. 2005> . the Sagittarius dSph 
JDinescu et al. 2005a> . the Fornax dSph (IPiatek et al. 20021 
IDinescu et al. 2004> . and the Carina dSph iPiatek et al. 2003> . 
These galaxies are all closer than 150 kpc and show motions 
between 0.2 and a few milliarcseconds (mas) per year. More 
distant galaxies, such as galaxies in the Andromeda subgroup 
at distances of ~ 800 kpc, have smaller angular motions, which 
are currently not measurable with optical telescopes. 

On the other hand, Brunth aler et al. (20 05) measured the 
proper motions of two groups of water masers on opposite 
sides of M33 at radio frequencies with the NRAO 1 Very Long 



Baseline Array (VLBA). A comparison of the relative proper 
motion between the two groups of masers and their expected 
motions from the known rotation curve and inclination of 
M33 led to a determination of a "rotational parallax" (730 
+ 168 kiloparsec) for this galaxy. This distance is consistent 
with recent Cepheid and tip of the red giant branch estimates 
JLee et al . 2002 McConnac hie et al. 20 05 ) and earlier distance 
estimates using the internal motions of water masers in the 
IC 133 region iGreenhill et al. 1993l|Argbn et al. 2004| i. 

Since the proper motion measurements were made relative 
to a distant extragalactic background source, the proper motion 
of M33 itself could also be determined. This measured proper 
motion of M33 is a first important step toward a kinematical 
model of the Local Group and was used to constrain the proper 
motion and dark matter content of the Andromeda Galaxy M3 1 
JLoeb et al. 2005> . 

Water masers in Local Group galaxies have also been found 
toward the Magellanic Clouds (e.g. IScalise & Braz 198TI and 
IC 10 (e.g. |Henkel, Wouterloot, & Bally 1986) . Other Local 
Group galaxies were searched, but no additional water masers 
have been detected (see Brunth aler et al. 20 06 and references 
therein). In this paper we report on VLBA observations of the 
maser in IC 10 to measure its motion. 

1.2. IC10 



Send offprint requests to: brunthal@mpifr-bonn.mpg.de The extragala cti c nature of IC 10 was first recognized by 

1 The National Radio Astronomy Observatory is operated by |Mayall (1935) 1 . |Hubble (1936| l proposed that it was likely 

Associated Universities, Inc., under a cooperative agreement with the a member of the Local Group and described it as "one 

National Science Foundation. of the most curious objects in the sky". However, ob- 
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servations of IC 10 have been always difficult because of 
the low Galactic latitude of 3°. IC 10 has been clas- 
sified as an IrlV galaxy (e.g. |van den Bergh 19 99 I, but 
IRicher et al. (2001) argue that it has more properties of a 
blue compact dwarf galaxy. It is also the nearest galaxy 
hosting a small starburst, evidenced by its large num- 
ber of Wolf-Rayet stars ( |Massey, Arm androff, & Conti 1992 1 
and the discovery of 144 HII regions JHodge& Lee 1990) . 
Observations of HI with the Westerbork Synthesis Radio 
Telescope by |Shostak & Skillman (1989) revealed that IC 10 
has a regularly rotating disk surrounded by a counter-rotating 
outer distribution of gas. 

The distance to IC 10 is subject to controversy because 
of difficulties caused by its low Galactic latitude. Early 
estimates claim a distance of 1-1.5 Mpc (Roberts 1962 1 
and 3 Mpc (Bottinelli, Gouguenheim, & Heidmann 1972 
Sandage & Tammann 1974 1. Huchtmeier (1979 > argued 
for a closer distance of 1 Mpc. The most recent deter- 
mination from multi-wavelength observations of Cepheid 
variables obtained a distance of 660+66 kpc to IC 10 
(Sakai, Madore, & Freedman 1999) , which we adopt through- 
out this paper. 

I C 10 hosts two know n H 2 masers, IC 10-SE and IC 10- 
NW (IBecker et al. 19931 . The strong SE-component was first 
detected by |Henkel, Wouterloot, & Bally (1986) and the 
whole spectrum of IC 10 showed strong variability since 
its discovery with flux densities between less than 1 Jy 
(Bec ker et al. 1 993 1 and a flare with a (single dish) flux density 
of 125 Jy ( IBaan & Haschick 1994i . Even intraday variability 
has been reported by [Argon et al. (1994) , but the strong com- 
ponent at vlsr ~ -324 km s~' has been persistent until now. 

2. Observations and Data Reduction 

We observed the usually brightest maser in IC 10-SE with the 
VLBA thirteen times between February 2001 and June 2005. 
The observations are grouped into six epochs, each comprising 
two closely spaced observations, except the first epoch with 
three observations, to enable assessment of overall accuracy 
and systematic errors (Table[Q. 

We observed in four 8 MHz bands in dual circular polariza- 
tion. The 128 spectral channels in each band yielded a channel 
spacing of 62.5 kHz, equivalent to 0.84 km s , and covered a 
velocity range of 107 km s _1 . The observations involved rapid 
switching between the phase-calibrator VCS1 J0027+5958 
from the VLBA Calibrator Survey ( Beasle yet al. 2002) , which 
is a compact background source with continuum emission, 
and the target sources IC 10 and NVSS J002108+591 132. 
NVSS J002108+591 132 is a radio continuum source from the 
NRAO VLA Sky Survey (NVSS) (ICondon et al. 19981 and is 
located only 8 arcminutes from the maser in IC 10. It was 
also detected in X-rays (Wang, Whitaker, & Williams 2005 1 
and is most likely also a background quasar. The redshifts 
of VCS1 J0027+5958 and NVSS J002108+591132 are not 
known. We switched sources every 30 seconds in the se- 
quence VCS1 J0027+5958 - IC 10 - VCS1 J0027+5958 - 
NVSS J002108+591132 - VCS1 J0027+5958 and achieved 
on-source integration times of ~ 22 seconds. The background 



Table 1. Details of the observations: Observing date, observa- 
tion length t b s , beam size 9 and position angle PA. 



Epoch 


Date 


tobslh] 


8 [mas] 


PA[°] 


I 


2001/02/09 


10 


0.53x0.33 


-15 


1 


2001/03/28 


10 


0.55x0.36 


-18 


I 


2001/04/12 


10 


0.63x0.37 


-5 


II 


2002/01/12 


10 


0.59x0.35 


-19 


II 


2002/01/17 


10 


0.64x0.32 


-22 


III 


2002/10/01 


10 


0.68x0.38 


-12 


III 


2002/10/11 


10 


0.61x0.34 


-5 


IV 


2003/12/12 


12 


0.52x0.33 


-15 


IV 


2004/01/10 


12 


0.50x0.33 


-23 


V 


2004/08/23 


12 


0.60x0.51 


-2 


V 


2004/09/18 


12 


0.54x0.35 


-17 


VI 


2005/06/01 


12 


0.60x0.50 


-11 


VI 


2005/06/07 


12 


0.56x0.39 


-6 



sources were assumed to be stationary on the sky. Since the 
phase-calibrator is separated by only 1° on the sky from the 
target sources, one can obtain precise angular separation mea- 
surements. 

From the second epoch on, we included geodetic-like ob- 
servations where we observed for 45 minutes 10-15 strong ra- 
dio sources (> 200 mJy) with accurate positions (< 1 mas) at 
different elevations to estimate an atmospheric zenith delay er- 
ror in the VLBA calibrator model (see Reid & Brunthaler 2 004 
and Brunthaler, Reid, & Falcke 2005 for a discussion). In the 
second and third epoch we used two blocks of these geodetic 
observations before and after the phase-referencing observa- 
tions. From the fourth epoch on, we included a third geodetic 
block in the middle of the observation. 

The data were edited and calibrated using standard tech- 
niques in the Astronomical Image Processing System (AIPS). 
A-priori amplitude calibration was applied using system tem- 
perature measurements and standard gain curves. Zenith de- 
lay corrections were performed based on the results of the 
geodetic-like observations. Data from the St. Croix station were 
flagged due to high phase noise in all observations. The maser 
in IC 10 and NVSS J002108+591 132 were imaged in AIPS. 
All detected maser features and NVSS J002108+591 132 were 
unresolved and fit by single elliptical Gaussian components. 

3. Results 

3.1. Spatial Structure 

In the first epoch, maser emission could be detected in 21 
channels spread over « 23 km s . The spatial distribution 
of the masers on 2001 March 28 can be seen in Fig. ^ 
It is similar to the distribution in earlier VLBI observa- 
tions of IC 10 by Arg on et al. (1994 |. The strongest compo- 
nent at a LSR velocity of « -324 km s _I is separated by 
« 10 mas or (projected) 6600 AU from the weaker com- 
ponents. This suggests that the emission is associated with 
a single object if the maser emission is similar to H2O 
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Fig. 1. Composite map of the H2O masers in IC 10 from 2001 
March 28. The area of each circle is proportional to the flux 
density of the respective component. The colors denote differ- 
ent LSR radial velocities with > -327 km s (red), -327 to 
-330 km s _1 (magenta), -330 to -338 km s _1 (green), and 
< -338 km s (blue). The crosses mark the positions of the 
maser emission detected by Argon et al. (1994i. The positions 
were aligned on the strongest maser component in the south- 
east. 



maser emission in Galactic star forming regions like W3(OH), 
W49 or SgrB2 (e.g. iReidet al. 19951 IWalkeretal. 19771 and 
|Kobayashi et al. 198 9 respectively). The weaker maser compo- 
nents form an apparent ring-like structure with a projected size 
of* 1.6 mas or 1060 AU. 



3.2. Variability 

The correlated flux density of the strong feature at -324 km s _1 
LSR velocity increased from 1.0 to 1 .5 Jy between the first two 
VLBA observations of the first epoch. The weaker components 
are also very variable and they can appear or disappear between 
observations. In the observations of the second epoch, the flux 
density of the strongest component was ~ 1.1 Jy and some of 
the weaker components disappeared. In the third epoch, we de- 
tected only the strong component with a flux density of ~0.7 Jy 
while the weak ring-like structure was not detected anymore. In 
the later epochs, the flux density of the remaining component 
dropped to ~ 0.2 Jy, 0.12 Jy, and 0.07 Jy in the fourth, fifth, and 
sixth epoch, respectively. 

3.3. Observed Motions 

The position offsets of the strongest maser feature in IC 10 are 
shown in Fig. [3] The uncertainties in the observations of the 
first epoch are larger than the others, because no geodetic-like 
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Fig. 2. The position of the maser in IC 10 relative to the phase- 
reference source VCS1 J0027+5958 in East-West (red trian- 
gles) and North-South (blue circles). The lines show a variance 
weighted linear fit to the data points. 
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Fig. 3. The position of NVSS J002 108 +591 132 relative to the 
phase-reference source VCS1 J0027+5958 in East-West (red 
triangles) and North-South (blue circles). The lines show a vari- 
ance weighted linear fit to the data points. 

observations were done to compensate the zenith delay errors. 
A rectilinear motion was fit to the data and yielded a value of 
-2+6 yuas yr~' toward the East and 20+6 juas yr~' toward the 
North. 

The position offsets of NVSS J002108+591 132 are shown 
in Fig. |3 A rectilinear motion was fit to the data and yielded a 
motion of -10+3 yuas yr _1 toward the East and -5+5 fias yr~' 
toward the North. Hence, NVSS J002108+591 132 shows a 
small but potentially significant motion in right ascension. The 
apparent motion of NVSS J002108+591132 could be caused 
by unknown systematic errors. 

The phase calibrator VCS1 J0027+5958 may have unre- 
solved structure, e.g. a core-jet structure. The observed image 
of the source is the convolution of the source structure and the 
synthesized beam of the VLBA. Flux density variations of the 
individual components could move the position of the observed 
image by a fraction of the beam size. Since the phase calibrator 
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Fig. 4. The position of the strongest maser feature in IC 10 rela- 
tive to NVSS J002108+591 132 in East- West (red triangles) and 
North-South (blue circles). The lines show a variance weighted 
linear fit to the data points. 



is assumed to be stationary, this would shift the positions of all 
target sources by the same amount. 

The observed motion of NVSS J002108+591132 could 
also be caused by some errors in the geometry of the corre- 
lator model (i.e. antenna positions, earth orientation param- 
eters). These errors would be similar for closely spaced ob- 
servations, but different for observations separated by sev- 
eral months. The angular separation between IC 10 and NVSS 
J002108+591 132 (8') is much smaller than the separation be- 
tween IC 10 and VCS1 J0027+5958 (1°), and the position shift 
induced by geometric errors would be similar for IC 10 and 
NVSS J002108+591 132. 

In both cases, the motion of the maser in IC 10 relative 
to NVSS J002108+591132 would be a better estimate of the 
proper motion of IC 10. However, it cannot be ruled out that 
the apparent motion of NVSS J002 1 08 +59 1 1 32 is caused by an 
unresolved core-jet structure in NVSS J002108+591 132 itself. 
Since strong sources are expected to show more jet-structure 
than weak sources, amplitude variations of VCS1 J0027+5958 
(70-290 mJy) are larger than those in NVSS J002108+591 132 
(6-11 mJy), and the angular separation between IC 10 and 
NVSS J002 108+59 1132 is much smaller than the separa- 
tion between IC 10 and VCS1 J0027+5958, we consider 
NVSS J002108+591132 as the better astrometric reference 
source. Fig.|4]shows the position of the strongest maser compo- 
nent in IC 10 relative to NVSS J002108+591 132. A rectilinear 
motion was fit to the data and yielded a motion of 6±5 /vas yr _I 
toward the East and 23+5 juas yr -1 toward the North and we 
will adopt these values for the proper motion of the maser. 

4. Discussion 

4.1 . Space Motion oflCW 

The measured proper motion v prop of the maser in IC 10 can 
be decomposed into a sum of several components, relative to a 
frame at rest at the center of the Milky Way: 



V prop — Vrot V pec + + Vic io (1) 

Here v ro , is the motion of the masers due to the internal 
galactic rotation in IC 10, v pec is the peculiar motion of the 
masers relative to circular galactic rotation and v is the appar- 
ent motion of IC 10 caused by the rotation of the Sun about the 
Galactic Center. The last contribution Vic lo is the true proper 
motion of the galaxy IC 10. 

The H2O masers in IC 10 are located within a massive HI 
cloud in the central disk. If one assumes that the masers are 
rotating with the disk, one can calculate its expected proper 
motion. Shostak & Skillman (1989 1 measure an inclination of 
45° from the ellipticity of its HI distribution. The masers are 
33 arcseconds (106 pc) east and 99 arcseconds (317 pc) south 
of the kinematic center. Unfortunately no position angle of the 
major axis was given. Wilcots & Miller (1998 1 used higher res- 
olution VLA observations of the HI content of IC 10 to fit a 
tilted ring model to the velocity field of the disk of IC 10. This 
model has a separate rotation speed, inclination and position 
angle for each ring. They find a highly inclined disk in the inner 
1 10 arcseconds with a position angle of « 75° and a rotational 
velocity of » 30 km s _1 . The position of the kinematic center 
of their tilted ring model was not given. If one combines the 
kinematic center of Shostak & Skillman (1989 1 with the incli- 
nation and position angle of Wilcots & Miller (1998 1, one gets 
an expected transverse motion for the maser (v rot ) of 26 and 1 1 
km s _1 toward the East and North, respectively. 

If one calculates the expected motion for different realis- 
tic scenarios (i.e. changing the kinematic center by +20 arc- 
seconds, and the inclination and the position angle of the ma- 
jor axis by +20°), one gets always transverse motions be- 
tween 20-30 km s _1 toward the East and 5-15 km s toward 
the North. The deviation of the motion of the masers from 
the galactic rotation is unknown. The radial velocity of the 
CO gas at the position of the maser in IC 10 is about -330 
km s -1 (Becker 1990 1, which is close to the radial velocity of 
the maser. In our Galaxy peculiar motions of star forming re- 
gions can be 20 km s as seen in W3(OH) (IXuet al. 2 006 
IHachisuka et al. 20061 . Hence, to be conservative, we adopt 
values of 25 ±20 and 10 ±20 km s~' toward the East and North, 
respectively. This translates to a ro , — 8 ± 6 and S ro , — 3 ± 6 
//as yr~' at a distance of 660±60 kpc. 

The rotation of the Sun about the Galactic Center causes 
an apparent motion of IC 10. The motion of the Sun can be 
decomposed into a circular motion of the local standard of 
rest (LSR) and the peculiar motion of the Sun. The peculiar 
motion of the Sun has been determined from Hipparcos data 
by |Dehnen & Binney ( 1998 1 to be in km s 1 U =10.00±0.36 
(radially inwards), Vo=5.25±0.62 (in the direction of Galactic 
rotation) and Wo=7.17±0.38 (vertically upwards). VLBI mea- 
surements of the proper motion of SgrA*, the compact ra- 
dio source at the Galactic Center, yield a motion of 6.379 
± 0.026 mas yr~' along the Galactic plane ( Reid et al. 1 999 
Reid & Brunthaler 20041 . Combined with a recent geometric 
distance estimate of the Galactic Center of 7.62 ± 0.32 kpc 
(Eisenhau er et al. 20 05 1, one gets a circular velocity of 225±10 
km s" 1 for the LSR. 
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This motion of the Sun causes an apparent proper motion 
of 38 ±4 fias yr~' in Galactic longitude and -6+1 ^as yr~' in 
Galactic latitude (for a distance of 660 kpc and Galactic co- 
ordinates of IC10 of / = 118.96°, b = -3.32°). Converted 
to equatorial coordinates, one gets a G = 37 + 4 ytzas yr~' and 
<S = -1 1 ± 1 //as yr . 

The true proper motion of IC 10 is then given by 

ff/c to = a prop - a rot - d 

= (6 (±5) - 8 (±6) - 37 (±4)) juas yr" 1 

= -39 + 9 juas yr 1 = -122 ± 31 km s" 1 
and (2) 

Sic io = 5 p r0 p - 5 m t - <5 

= (23 (±5) - 3 (+6) + 1 1 (±1)) //as yr" 1 

= 31 ± 8 yuas yr 1 = 97 + 27 km s" 1 

The measured systematic heliocentric velocity of IC 10 
(-344 + 3 km s^. lde Vaucouleurs et al. 199H is the sum of the 
radial motion of IC 10 toward the Sun and the component of the 
solar motion about the Galactic Center toward IC 10 which is 
-196+10 km s -1 . Hence IC 10 is moving with 148+10 km s" 1 
toward the Sun. 

The proper motion and the radial velocity combined give 
the three-dimensional space velocity of IC 10. The total veloc- 
ity is 215+42 km s _1 relative to the Milky Way. This velocity 
vector is shown in the schematic view of the Local Group in 
Fig. 13 Here, we used Cartesian coordinates, where the Sun 
is located at the origin and the Galactic Center is located at 
(x,y,z)=(7.62,0,0) (see AppendixlAlfor details). 

4.2. Local Group Dynamics and Mass of M31 

If IC 10 or M33 are bound to M31, then the velocity of the 
two galaxies relative to M3 1 must be smaller than the escape 
velocity and one can deduce a lower limit on the mass of M3 1 : 

M M 3i > (3) 

A relative velocity of 147 km s - for a zero tangential 
motion of M31 - and a distance of 262 kpc between IC 10 
and M31, gives a lower limit of 6.6 xl0 n M Q . One can repeat 
this calculation for any tangential motion of M31. The results 
are shown in Fig.|6](top). The lowest value of 0.7 xlO n M is 
found for a tangential motion of M31 of -130 km s toward 
the East and 35 km s _I toward the North. 

For a relative motion of 230 km s _1 between M33 and 
M31 - again for a zero tangential motion of M31 - and a 
distance of 202 kpc, one gets a lower limit of 1.2 xlO 12 M 
( Brunthal er et al. 20 05 1. Fig.[6](top) shows also the lower limit 
of the mass of M3 1 for different tangential motions of M3 1 if 
M33 is bound to M31. The lowest value is 4 xlO n M for a 
tangential motion of M31 of -1 15 km s _1 toward the East and 
160 km s" 1 toward the North. 

|Loeb et al. (20 05 1 find that proper motions of M3 1 in nega- 
tive right ascension and positive declination would have lead to 




Fig. 5. Schematic view of the Local Group from two viewing 
angles with the space velocity of IC 10 and M33 (the latter 
taken from Brunth aler et al. 2005> and the radial velocity of 
Andromeda relative to the Milky Way. The blue cross marks 
the position of the Local Group Barycenter (LG BC) according 
to |van den Bergh (19~99) . 

close interactions between M31 and M33 in the past. These 
proper motions of M31 can be ruled out, since the stellar 
disk of M33 does not show any signs of strong interactions. 
|Loeb et al. (2005) used a total mass of M31 of 3.4xlO 12 M in 
their simulations. Although simulations with lower masses of 
M31 yield weaker interactions, motions in negative right as- 
cension and positive declination are still ruled out. 

Thus, we can rule out these regions in Fig. [6] This results 
in a lower limit of 7.5x10' 'M for M31 and agrees with a 
recent estimate of 12.3*g 8 x 10 11 M derived from the three- 
dimensional positions and radial velocities of its satellite galax- 
ies JEvans & Wilkinson 2000V 

5. Summary 

We have presented astrometric VLBA observations of the H2O 
maser in the Local Group galaxy IC 10. We detected a ring- 
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Fig. 6. Top: Lower limit on the mass of M3 1 for different tan- 
gential motions of M31 assuming that M33 (dashed) or IC 10 
(solid) are bound to M31. The lower limits are (4, 5, 7.5, 10, 
15, 25)xlO n M G for M33, and (0.7, 1, 2.5, 5, 7.5, 10, 15, 
25)xlO n M for IC 10, rising from inside. The colour scale in- 
dicates the maximum of both values. Bottom: The colour scale 
is the same as above and gives the lower limit on the mass of 
M3 1 . The contours show ranges of proper motions that would 
have lead to a large ammount of stars stripped from the disk of 
M33 through interactions with M31 or the Milky Way in the 
past. The contours delineate 20% and 50% of the total number 
of stars stripped (Loeb et al. 2005 1. These regions can be ex- 
cluded, since the stellar disk of M33 shows no signs of such 
interactions. 

like structure in one epoch with a projected diameter of ~1060 
AU. We measured the proper motion of the maser relative to 
two background quasars. Correcting for the internal rotation 
of IC 10 and the rotation of the Milky Way this measurement 
yields a proper motion of -39+9 yuas yr~' toward the East and 
31+8 fj. as yr~' toward the North, which corresponds to a total 
space velocity of 215+42 km s _1 for IC 10 relative to the Milky 
Way. If IC 10 and M33 are bound to M31, one can calculate a 
lower limit of the mass ofM31 of 7.5 xlO n M . 

Acknowledgements. This research was supported by the DFG Priority 
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Appendix A: Coordinate Transformations 

We define a Cartesian coordinate system where the Sun is lo- 
cated at the origin. To convert from Galactic to Cartesian coor- 
dinates, we used 
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R sin I cos b 



(A.l) 



where R is the distance, I is the Galactic Longitude, and b 
is the Galactic Latitude. 

To convert the proper motions in right ascension (a) and 
declination (6) into proper motions in Galactic Longitude (/) 
and Galactic Latitude (b), we used 
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where 6 is the required rotation angle at the position of the 
source (i.e. 6, cw = -6.94°, 6 M33 = -348.95°, M 3i = -2.74°). 

To compute the velocities in Cartesian coordinates, we con- 
structed for each source three orthogonal vectors e r , e\, and ei, 
given by 
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(A3) 



These vectors can be normalized to give three orthogonal 
unit vectors e r , e), and e~b- Then the total velocity v tot of an ob- 
ject with proper motions /, b, and the radial velocity v ra d is 



Vtot = 1 ei + b e h + v rdd e r . 



(A.4) 



